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Abstract We propose that relativistic Galactic jets like those observed in GRS 1915+105 and GRO J1655-40 
may produce a small but measurable contribution to the cosmic ray (CR) spectrum. If these jets contain cold 
protons and heavy ions (as in the case of SS433), it is plausible that this component will consist of a narrow 
spectral feature, with a mean particle energy corresponding roughly to the bulk kinetic particle energy in the 
beam, rj e tm p c 2 . Based on the current estimates of Tjct, this feature will fall into the range of 3-10 GeV. The 
presence of several sources with different Tjet will lead to the superposition of several such peaks. In addition to 
the narrow peaks, diffusive particle acceleration should also produce a powerlaw, whose low energy cutoff at or 
above Ff et m p c 2 would be visible as an additional spectral feature. The large metallicities measured in several 
binary companions of jet sources suggest that this CR component could have an anomalous composition compared 
to the bulk Galactic CR spectrum. We provide estimates of the effects of adiabatic losses which are the greatest 
challenge to models of narrow band CR production in microquasar jets. While the total energy contained in the 
microquasar CR component is highly uncertain, the local CR spectrum in the vicinity of any microquasar should 
be severely affected. The upcoming AMS 02 experiment will be able probe the low energy CR spectrum for such 
components and for composition anomalies. The spectrally peculiar gamma ray emission produced by interaction 
of the ISM with CRs surrounding microquasars might be detectable by GLAST. If the presence of a microquasar 
CR proton component can be ruled out observationally, this argument could be turned around in favor of electron- 
positron jets. We show that existing OSSE/GRO and future INTEGRAL data on the Galactic 511 keV line flux 
put interesting constraints on the particle content of microquasar jets. The process of CR production in relativistic 
flows inside the Galaxy is fundamentally different from the standard picture of CR production in nonrelativistic 
shocks in supernova remnants, because the particles injected by a relativistic flow are already relativistic, without 
any need for diffusive acceleration. 

Key words, acceleration of particles - ISM: cosmic rays - ISM: jets and outflows - shock waves - black hole physics 
- Gamma rays: theory 



1. Introduction 

During the past decade, a new class of relativistic jet 



jets most likely operate at a variety of Lorentz factors, 
which are currently not well determined. These ejections 
can contain energies in excess of 10 44 ergs ( Mirabel & 



source ^ has been discovered in Lhc Galaxy: GalacLicsuper- Rodriguez Tggj [Hjcllming fc R.upcnT99l [Fender et al 



lumin a^ radio sotnees (ot niicrociuasars), exeinpimed by mg ^ and ^ occur at a Galactic rate of f 
the pr ototypical sources GRS 1915+105 and GRO J1655- ^ 

40. Radio outbursts from these transient X-ray binaries, 
containing accreting black holes or neutron stars, are asso- 
ciated with ejections of plasma at relativistic bulk veloci- 
ties. The bulk Lorentz factors of these flows inferred from 



observations are of the order of rj et ~ 3 — 10, though exact 
observational determination of Tj Ct is difficult and in any 
case it is unlikely that a universal value of Tj Ct holds for all 
jets. For reasons of simplicity and for lack of better knowl- 
edge we will adopt a fiducial value of Tj Ct ~ 5 for numerical 
examples, with the explicit understanding that Galactic 



A large part of the kinetic energy transported by these 
jets is transferred into random, isotropic particle energy at 
the interface between the jet and the ambient medium, the 
working surface. Because the jets are relativistic, the parti- 
cles leaving the working surface must a priori be relativis- 
tic without any need for diffusice acceleration. This mech- 
anism of accelerating relativistic cosmic ray ( CR ) particles 
is fundamentally different from CR production in the non- 
relativistic shocks of supernova remnants (SNRs), which 
provides the bulk of the Galactic CRs. 
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While the momentum gain for particles crossing a non- 
relativistic shock is small (of order Sp/p ~ v/c), the large 
momentum gains encountered in relativistic shocks (of or- 
der Sp/p ~ r, where T is the shock Lorentz factor) should 
lead to the formation of distinct spectral features in the 



spectrum (see Sect. 3.3.1). Thus, unlike the CRs produced 
in SNRs, which follow a smooth powerlaw spectrum, the 
CRs produced in relativistic flows, like those encountered 
in microquasars, should show clearly distinguishable, and 
possibly narrow, spectral features. 

If these particles can escape the working surface with- 
out suffering significant adiabatic energy losses, they will 
diffuse through interstellar space, and will thus contribute 



to the Galactic cosmic ray (CR) spectrum 



Based on the arguments presented in this paper, we 
conclude that an additional component of CRs generated 
by relativistic jets in microquasars should exist in the 
Galaxy. Initially, this component should consist of nar- 
row peaks, with peak energies corresponding to rj ct m p c 2 
from different jet sources. 

There are many mechanisms which might broaden 
these features. However, any observational limits on their 
existence would give us additional information about the 
physics of microquasar jets and the physical conditions in 
relativistic shocks. Below, we will discuss the main mech- 
anisms which could smooth out the component under dis- 
cussion. In this paper, we content ourselves with present- 
ing order of magnitude estimates only, since the goal of the 
paper is to point out to the CR community that, in ad- 
dition to CR acceleration in supernova remnants (SNRs), 
there is another very effective mechanism to release rela- 
tivistic particles in the Galaxy. Traces of these particles 
might be hidden in the observed CR spectra, in 7-rays 
with energies of a few 100 MeV, and possibly in electron- 
positron annihilation line emission from regions close to 
the location of microquasars in the Galactic plane. 



2. Outline of the model 

In this section we will present a general description of the 
model proposed for CR production in microquasars. 

As we will argue below, it is likely that the plasma 
traveling far downstream in the jet towards the working 
surface is cold (the mean particle velocity is (v 2 ) <C c 2 in 
the rest frame of the jet plasma), especially if microquasar 
jets are composed of electron-ion plasma. The same is, of 
course, also true for the undisturbed interstellar medium 



1 The production of ultra-high energy CRs in Galactic X-ray 
binar ies has been discussed in the past in the context of Cyg 



X-3 (Lloyd-Evans et al. 1983; 



Manchanda et al. 1993 



Mitrej 



1994), though recent studies have not confirmed the pr esence 



of ultra-high energy gamma rays ( Borione et al. 1997 ). This 
scenario is shown as a cartoon in Fig. |l| The production of 
CRs in microquasars has so far only received minor attention in 
the literature (for example, it w as mentioned as an alter native 
source to produce all CRs by 



(ISM). Thus, the bulk of the plasma transported to the in- 
terface between the jet and the ISM (for simplicity we will 
call this interface the working surface of the jet, regardless 
of its detailed physical structure) is initially cold. 

This conjecture is inspired by observations of the 
mildly relativistic jets in SS433 (the best studied relativis- 
tic Galactic jet to date, albeit mildly relativistic and not 
considered a microquasar). In this source, red- and blue- 
shifted Balmer H/3 and other optical recombination lines, 
usually radiated by plasmas with temperatures of order 
10 4 K ~ 1 eV, allow the determination of the bulk veloc- 
ity of the flow: 0.26 c. This velocity is remarkably constant 
over the 20 years the source has been observed flMargon 
1984 |Mifgrom et al. 198%). ASCA (|Kotani et al. 1998j ) and 



requires further study. 



Dar & De Rujula 2001), and 



recent Chandra ( Marshall et al. 2002| ) observations of X- 
ray lines of hydrogen- and helium-like ions of Iron, Argon, 
Sulfur, and Oxygen show that these ions are moving in 
the flow with the same velocity, 0.26 c. This X-ray emitting 
plasma at temperatures of T ~ 10 7 - 10 s K - 10 3 - 10 4 eV 
is observed at much smaller distances (~ 10 11 cm) from 
the central compact object than the optical line emission 
region (~ 10 14 cm). A striking feature of the SS433 jet is 
that the plasma is observed to be moving with relativis- 
tic velocities. Yet, at the same time the jet plasma itself 
shows very little line broadening (i.e., it is cold). 

If microquasar jets are similar to the SS433 jets in com- 
position and properties (i.e., cold electron- ion plasma at 
relativistic bulk speeds) the consequences for the interpre- 
tation of these jets will be far reaching, as we will argue 
below. Independent from this argument, the radio syn- 
chrotron emission detected from microquasar jets and sev- 
eral radio nebulae surrounding microquasar sources (see 
Sect. 3.1) is clear evidence for the presence of relativistic 
electrons, which, when released into the ISM, will act as 
cosmic ray electrons. 



2.1. The standard picture for jet working surfaces 

The standard picture for the interface between powerful 
radio galaxies and their environment is a strong double 
shock structure (forward shock into the ISM and reverse 
shock into the jet), shown in Fig. ||. The shocked jet mate- 
rial is shed at the head of the jet and inflates an enshroud- 
ing cocoon around the jet, filled with relativistic plasma, 
which has gone through the terminal shock. Such a sce- 
nario might also be relevant for the terminus of Galactic 
relativistic jets. A similar picture arises if the jets are com- 
posed of discrete ejections, propagating into an external 
medium at relativistic speeds, as sketched in Fig. |3|. 

A cold upstream particle crossing an ultra-relativistic 
shock into a downstream region with relative Lorentz fac- 
tor r ro i ~ Tjot will have an internal energy of 7771c 2 ~ 
Tjetmc 2 in the downstream frame after the first shock 
crossing. Consequently, all initially cold particles will leave 
the shock with about the same specific energy Fj Ct c 2 . 
Particles can pick up additional energy if they cross the 
shock multiple times, which is the basis of diffusive shock 
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Figure 2. Cartoon of the standard picture of the interface 
between jet and ISM (left panel), as envisaged to apply in 
FR II radio galaxies. The injection of relativistic particles 
can occur either in the reverse or forward shock. Right 
panel: cartoon of particle trajectories for particles cross- 
ing the shock only once (upper solid line) and particles 
participating in diffusive shock acceleration (lower solid 
line), particle scattering indicated as stars. 



Figure 1. Cailuun uf the proposed mudul uf GR piuduc- 



tion 



parti cle energy is just Fj Ct mc 2 (e.g. Blandford fc McKee 



m miciuquasais. llic mleilaci! uulwouii liie lulalivis- 

tic jet and the ISM is a natural site for the production 
and release of relativistic particles. 



acceleration schemes like Fermi acceleration, resulting in 
the formation of a powerlaw distribution. However, as has 
recently been shown by Achterberg et al. (2001 ) , the bulk 
of the particles crossing a relativistic shock escape after 
the very first shock passage and will therefore not partic- 
ipate in diffusive shock acceleration. It is these particles 
that carry off the bulk of the dissipated jet energy. 

As a result, the bulk of the particles might leave the 
shock with a narrow energy distribution, peaking at an 
energy close to the specific kinetic energy of the jet: 
(7mc 2 ) ~ Tjet mc 2 , with an energy width similar to or 
higher than the Lorentz transformed thermal velocity, 
A7 ~ 2 7c s (i.e., very narrow, since the internal sound 
speed c s is small: Cs < 1). 

Whether this narrow distribution will be preserved as 
the particles travel away from the shock, or whether it 
will be thermalized, depends on the efficiency of collective 
plasma effects and small angle scattering on magnetic field 
irregularities, which are also needed to isotropize the parti- 
cle distribution. If collective effects are strong, the particle 
spectrum will be broadened into a relativistic Maxwell- 
Boltzmann distribution, with a temperature correspond- 
ing to the value given by the relativistic Rankine-Hugoniot 
jump conditions. In the case of a strong, ultra-relativistic 



1976). In this case the relativistic proton plasma in the 



shocked ISM is equivalent to the X-ray emitting gas in 
SNR shocks. However, in microquasars shocks we have 
extremely rarefied, relativistic particles with a relatively 
narrow thermal (i.e., not powerlaw) energy distribution. 

2.2. Termination without strong shocks? 

However, the structure of relativistic shocks is still not 
well understood and it might be that this interface is not 
a simple double shock structure. It could be significantly 
different in nature. For example, the jet could be mag- 
netically connected with the environment, i.e., if the flux 
tubes join smoothly with the large scale magnetic field of 
the ISM, as shown in the car toon in Fig. [| (note, how- 
ever, that Lubow et al. 1994 showed that realizing such 



shock, this is simply kT ~ l/3Tj ct mc , i.e., the mean 



configuration is rather difficult). 

In such a case the shock would be replaced by stochas- 
tic pitch angle scattering of the particle distribution (this 
can occur if the jet is moving sub-Alfvenically, for exam- 
ple). Since the jet plasma is traveling relative to the ISM, 
such a scenario would excite strong two-stream instabili- 
ties at the interface between ISM and jet plasma, which 
would isotropize and possibly thermalize the particle dis- 
tribution of the jet very quickly. The deposition of jet 
thrust would then imply that this interface is itself mov- 
ing through space. Precession, as observed in SS433 (e.g., 
Milgrom 1979|) and suggested to be present in GRO 1655- 
40 (Hjcllming & Rupcn 1995), will significantly alter the 
dynamical balance between ISM and jet plasma, as will 
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Figure 3. Left: Cartoon of a jet composed of discrete ejec- 
tions with precession. In such a non-stationary picture, 
each ejection is slowed by its interaction with the ISM 
(which might be disturbed by previous ejections). This in- 
teraction will likely happen in the form of a forward shock 
(into the ISM). ISM particles will leave the shock with en- 
ergies of order rj ot m p c 2 (here, Tj ot is the Lorentz factor of 
the individual ejection, which decreases with time). Right: 
Cartoon of particles crossing a relativistic forward shock, 
likely the appropriate scenario for a cold ejection running 
into the ISM. Same nomenclature as in Fig. 0. 























ISM 





Figure 4. Cartoon of particle injection by a jet without a 
strong shock at its interface with the ISM (e.g., if the mag- 
netic structure of the jet is connected with the ISM and 
the jet is sub-Alfvenic). Dots represent CRs transported 
by the jet and released into the ISM at the bulk speed 
of the jet, grey arrows represent magnetic flux tubes. The 
right flux tube representation depicts the case where adia- 
batic losses are important, the left shows a case where the 
field flux tube is smoothly fused with the stochastic inter- 
stellar field and adiabatic losses might not be dominant in 
particle transport. 



the time dependent nature of the interface if the jets are 
composed of discrete ejections. 

If furthermore the magnetic field is stochastically 
tangled on small scales, the detailed behavior of the 
plasma could be very complicated, with a gradual change 
from relativistic, ballistic motion to random propagation. 
Qualitatively, this would be comparable to extragalactic 
FR I sources (though the exact nature of the dynamics in 
FR I sources is not yet clear, either). 

In such a case, the absence of a strong shock would 
preclude diffusive shock acceleration (though stochastic 
acceleration might still exist if particles scatter off of rel- 
ativistic turbulence which might exist in the transition 
region between jet and ISM). Only the narrow or thermal- 
ized component with mean energy of Tjot m p c 2 and strong 
cutoff at higher energies would exist. 

2.3. Spectral characteristics 



Based on Sect. 2.1 and Sect. 2.2, we therefore predict that 



each Galactic microquasar produces a narrow component 
of CRs, which peaks at an energy T- iet mc 2 . For protons, 
this energy should fall into the range of rj et m p c 2 ^3 — 10 
GeV. The superposition of these spectral signatures from 
several microquasars will appear as a broad feature in the 
energy range from 3 — 10 GeV. 

Even if most of the particles are thcrmalized down- 
stream, the spectrum will still show a steep turnover be- 
yond energies of order 3kT ~ rj ct m p c 2 (see the dashed 
curve in Fig. |^), which will appear as an edge-like feature 
in the overall CR spectrum. 



Similarly, a number of other processes will tend to 
broaden any narrow component produced in the working 
surface, including adiabatic losses (competing with diffu- 
sion of particles out of the loss region, see Appendix ^ 
and the insert in Fig. ||) and solar modulation. The effect 
of these processes will be to spread particles to lower en- 
ergies, leaving the strong turnover/cutoff above energies 
of Tjet m c 2 intact. 

The only serious cooling these CR protons at energies 
of a few GeV might experience will be adiabatic losses, 
which will occur if the particles are confined to an ex- 
panding plasma volume (e.g., if it is overpressured with 
respect to the environment). However, since many pro- 
cesses can lead to increased diffusion of these particles, it 
appears plausible that a large fraction of the CRs might 
escape before they suffer strong adiabatic losses. 

If a component of cold electrons is also present in 
the jets in addition to the observed powerlaw electrons, 
a similar, very low energy relativistic electron component 
(around 2-5 MeV) might appear. However, it would con- 
tain only a fraction m e /m v of the energy in the proton 
component. 

The remaining fraction of particles (both ions and elec- 
trons) which do not escape the shock after the first shock 
passage and thus perform multiple shock crossings will 
be accelerated diffusively to a powerlaw-like distribution. 
Only the high energy tail of this powerlaw-like electron 
component is directly observable via synchrotron radio 
emission. 

Likening the acceleration of particles crossing a rela- 
tivistic shock to the problem of Compton up-scattering of 
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Figure 5. Left Panel: Sketch of the predicted contribution from a microquasar to the Galactic CR spectrum for 
Tjct = 5. a) Dotted line: narrow component, corresponding to cold (T ~ 7 x 10 10 K) upstream particles simply 
isotropized upon crossing the working surface (Lorentz transformed energy distribution), arbitrary normalization; b) 
dashed line: Maxwell-Boltzmann component from particles thermalized in the shock (normalized to contain the same 
power as component a); c) dash-dotted line: powerlaw-like component (escape probability 99%, normalized to contain 
10% of component a) for the case of strong scattering. Note the similarity to Comptonization of low energy seed 
photons by a relativistic medium with low optical depth (i.e., high escape probability); d) dash-triple- dotted line: 
powerlaw spectrum for weak scattering case (Achtcrberg ct al. 2001) for the same normalization as component c). 
Note the difference in the position of the first peak between c) and d), ~ r? et vs. ~ r? ot . Right Panel: estimate of the 
effects of adiabatic losses on the particle distribution for different ratios of the adiabatic losses time T ac j of the particles 
to the escape time r osc , compared to the injected particle distribution (solid) line for a pressure differential between 
shock and ISM of p s hock/piSM = 3 x 10 8 (see Appendix [bJ) : strong adiabatic losses (T a d/r csc = 10~ 3 , dashed line), 
intermediated losses (r a d = t csc , dash-dotted line), and weak adiabatic losses (T a d/T eS c = 10 3 , dotted line). 



low energy photons on relativistic thermal electrons (see, 
for example, Pozdnyakov et al. 1983), we note that a par- 
ticle scattered both up-stream and down-stream of the 
shock will experience an energy gain by a factor of order 
^rei P er cross i n g cycle, as was argued by, where r re i ~ Tj Ct 
is the relative Lorentz factor betwee n upstream a nd down- 
stream plasma. This was argued by Victri (1995 ), applied 
to the acceleration of particles in gamma ray burst shocks. 
This will lead to the production of several peaks in the 
spectrum. The input spectrum for this up-scattering pro- 
cess is the narrow particle population produced in the ini- 
tial shock crossing (discussed above), and thus peaks will 



appear at energies ~ few x rj 0t 



2i+l 



TO D C 



where i is the 



number of shock crossing cycles performed by the particle. 
The normalization of each peak, and thus the approximate 
powerlaw index, is determined by the escape probability 
of the particles (similar to the optical depth in inverse 
Compton scattering). The resulting spectrum is sketched 
in Fig. |. 



Note, however, that Achterberg et al. (2001) argue that 
higher order shock crossings do not lead to energy gains of 
order rjr ol . In their treatment, scattering is limited to very 
small angles and the energy gain is only of order unity, 
and thus the position of the peaks would be much more 
closely spaced, resembling a powerlaw much more than in 
the Compton scattering analogy discussed in the previous 
paragraph. The low energy turnover (or cutoff) of this 



powerlaw distribution would then be located roughly at 
r? et ^ 10— 100 GeV. At higher energies, multiple scatter- 
ing will form a powerlaw with indexQ s ~ 2.3. According 
to this simple approach, the difference between these two 
pictures is therefore the energy of the second peak (~ r? et 

vs. ~ r? t ). 

Since the structure of relativistic shocks, and their 
presence in the working surfaces of microquasar jets are 
subject to considerable uncertainty, the observational dis- 
covery of any of the features discussed in this paper (and 
in particular the second peak, which would help to dis- 
tinguish between the two scenarios of diffusive accelera- 
tion mentioned in the previous paragraphs, see Fig. ||) or 
evidence of their absence would be important input into 
theories of relativistic shocks. 



2.4. Comparison with the canonical CR powerlaw 
component 

The CR components described above and shown in Fig. || 
should be compared to the well known powerlaw CR 
component observed near earth: CR protons and nuclei 
have a powerlaw spectrum with a uniform slope around 



Note that powerlaw indices from cosmic ray modified 
shock s tend to be steeper at l e ast near the low energy cutoff 



(e.g., Berezhko fc Ellison 1999; Ellison fc Double 2002) 
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s ~ 2.5 — 2.7 over an extremely broad energy range from 
1 GeV up to 10 15 eV. Locally, the CR energy density is of 
order 10~ 12 ergs cm -3 . Observations of the light elements 
produced by spallation reactions show that the lifetime 
t cr, disk of relativistic protons in the Galactic disk is of or- 
der rca, disk ~ 1-5 x 10 7 yrs ( Yanasak et al. 2001), while the 



lifetime in the Galactic halo tcr, halo is close to 10 years 



(e.g., Ginzburg 1996). This enables us to estimate the CR 
luminosity of the Galaxy as Lcr ~ 4x 10 40 ergs s _1 , with 
a diffusion coefficient n inside the disk of 



that there is an unknown acceleration mechanism with 
100% efficiency of transforming of the mechanical beam 
energy into a relativistic powerlaw distribution. Given 
these premises, we can state that one of the following two 
statements must hold: 1 ) Either an additional hadronic 
CR component exists (though it may be so weak that de- 
tection inside the solar system is impossible) or 2) All jets 
are electron-positron dominated (in which case an addi- 
tional CR electron-positron component should exist). 



H, 



disk 2 



TCR, disk 



2 X 10 



2<S ' 



kpc 



10 



28.3 



TIB 



k-28. 3 (1) 2.5. Peculiar abundances 



where -f/kpc is the disk height in kpc, rig is the CR lifetime 
in the disk in units of 15 Myrs, and «28.3 is the ISM diffu- 
sion coefficient, normalized to a value of 2 x 10 28 cm 2 s _1 . 
The current paradigm for the bulk of the CRs observed 
in the vicinity of the earth is that they are produced by 
shock acceleration in the decelerating blast waves of SNRs 
(e.g. Blandford & Ostriker 1978). The usual assumption 
is that about 5% of the mechanical energy of the SNR are 
converted into CR energy. 

There is no doubt that in the vicinity of an active 
microquasar the low energy part of the Galactic CR spec- 
trum must be strongly distorted. As a result, smooth max- 
ima or edge-like features should exist in the few GeV 
range of the CR spectrum. For a distant observer, the sig- 
nals from several sources will be superimposed due to the 
long diffusion time through the galaxy. Integrally, though, 
deviations from the powerlaw spectrum expected in diffu- 
sive shock acceleration models should be observable. 

Energy estimates which we present below show that 
this CR component produced in microquasars might con- 
tribute measurably to the spectrum of the CR protons 
in the energy band mentioned above. We will argue that, 
globally, microquasars should contribute upward of 0.1% 
of the total Galactic CR power. However, the locally mea- 
sured (i.e., near earth) relative strength of the proposed 
CR components produced in microquasars compared to 
the canonical CR powerlaw distribution is highly uncer- 
tain, as it depends on the history of microquasar activity 
in our Galactic neighborhood. 

Given these uncertainties, it might be rather difficult 
to detect the tiny deviations in the CR spectrum caused by 
distant microquasars (further complicated by the strong 
effects of solar modulation at and below the predicted en- 
ergy range). However, they might be measurable by the 
upcoming AMS 02 experiment (e.g. Barrau 2001 ), which 
will offer unprecedented sensitivity and will be launched 
during the upcoming solar minimum (reducing the effects 
of solar modulation significantly). Traces of such a compo- 
nent might also be present in already existing high quality 
data sets from past or ongoing experiments, such as IMAX 



( pVIenn et al. 2000[ ) or CAPRICE ( [Boczio ct al. 1999| ). 

Absence of any traces of spectral deviations in the up- 
coming AMS 02 experiment might become a strong argu- 
ment in favor of electron-positron jets in Galactic superlu- 
minal radio sources or, alternatively, it would demonstrate 



Recent studies of the chemical composition of Galactic 
transients indicate that the binary companions of several 
microquasar sources are metal enriched: the heavy ele- 
ment abundances (e.g., N, O, Ca, Mg) in the optical coun- 



terparts of the X-ray sources GRO J 1655-40 (Israeliari 
|ct al. 1999[ ) and V4641 Sgr ( prosz et al. 2001 ) exceed 



solar abundances by about an order of magnitude, much 
more so than t he observed ove rabundance of the same el- 
ements in CRs (IZombcck 1990D - Furthermore, the relative 
abundances between these elements are unusual compared 
to either solar abundances or the abundance pattern ob- 
served in the bulk of the CR spectrum. 

These abundance anomalies in GRO J1655-40 and 
V4641 Sgr could be the result of mass exchange between 
two rapidly evolving massive stars or enrichment of the 
normal stellar atmosphere during the supernova explosion 
of the primary predecessor. Accretion brings these abun- 
dance anomalies into the jet creation region in the inner 
disk, from where they could be transported out by the 
jet, eventually producing CRs by the mechanism outlined 
above. Similarly, Cyg X-3 is known to have an extremely 



hydrogen deficient Wolf-Rayet companion (van Kerkwijk 
|ct al. 1992| , |1996| ; Render et al. 1999b| ), which could also 
lead to a large overabundance in helium and heavier ele- 
ments relative to hydrogen in the produced CR spectrum. 

Therefore, Galactic jets might be responsible for part 
of the observed CR abundance anomalies. Moreover, the 
CR component produced in relativistic jet sources inside 
the Galaxy might show rather unusual chemical abun- 
dances, in comparison with the bulk of the CRs in the 
powerlaw population. This would immediately distinguish 
Galactic jets from other CR creation mechanisms. The 
comparison between the measured abundances in the en- 
ergy range where we expect Galactic jet sources to con- 
tribute (of order a few GeV) with those measured in the 
pure powerlaw regime will thus be an important probe to 
search for the proposed CR component. Note that the CRs 
produced in SNRs originate in the external shock of the 
swept up ISM, thus the abundances of the produced CR 
spectrum reflect the ISM, which might have been enriched 
by a pre-collapse wind, but will not show the peculiar 
abundance of the SN ejecta. Because all the spectral com- 
ponents accelerated in microquasars originate from the 
same plasma, they should all show the same abundance 
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pattern. This could be a way to associate spectral features 
at different CR energies with a microquasar origin. 

A second way to distinguish particles accelerated in 
the relativistic shocks of Galactic microquasars from those 
accelerated in non-relativistic SNR shocks is the different 
energy — particle mass relation: All particles in relativistic 
cold jets have the same Lorentz factors. Since single-pass 
shock acceleration will accelerate all particles to roughly 
the same random Lorentz factor, the peak energy for dif- 
ferent species will be proportional to their rest mass (i.e., 
a fixed energy per nucleon) . Electromagnetic acceleration 
processes would instead produce particle energies propor- 
tional to Z/A. This difference might again be measurable 
by AMS 02, and might already be present in CR data on 
heavy nuclei from experiments like HEAO-3 flEngelmann 
et al. l|990|) or ACE (pinns et al. 200l|). 



power in excess of L^ n ~ few x 10 — few x 10 ergs s 



3. Energetics and spectra of microquasar CRs 

3.1. Microquasar energetics 

The best studied cases of microquasar activity show 
that large amounts of kinetic energy can be liberated: 
Conservative equipartition estimates of the energy re- 



( [Marshall et al. 2002| ; |Margon 1984 ppencer 1984| ), which 
is already of order 1 - 10% of the total Galactic CR lumi- 
nosity. While the jets in SS433 are only mildly relativis- 
tic, and the production of observable, relativistic CRs thus 
falls under similar restrictions with regard to particle ac- 
celeration efficiency as supernovae, the example of SS433 
does show that Galactic jets are capable of releasing im- 
pressive amounts of kinetic energy. Thus, jets from objects 
like SS433 (with mildly relativistic bulk speeds) might be 
an important source of sub-CRs in the Galaxy, influencing 
heating and ionization of the ISM. 

Since the subject of Galactic microquasars is still 
relatively young, and many of the known sources have 
only been discovered in recent years, estimating the true 
Galactic rate of radio outburst events and thus the to- 
tal Galactic power in relativistic jets is difficult. Taking 
the interval from 1994 through 2000, there were at least 7 
well observed giant radio outbursts comparable in strength 
with GRS 1915+105 ( corrected for Galactic dista nce) in 
the source s Cyg X-3 (|Mioduszewski et al. 2001] ), GRO 
J1655-40 (|Hjcllming fc Rupen 1995[), GRS 19 15+105 
([Mirabel fc Rodriguez 1994|; [Fender et al. 1999a|), V4641 



leased in the major outbursts of GRS 1915+105 give 



E kin > 2 x 10 44 ergs QMirabel fc Rodriguez 1999|; [Fender 



Sgr (Orosz et al. 2001), and XTE J 1748-288 (Fender fc 



et al. 1999a), released over a period of a few days at most. 

Existing radio m onitoring data ( Pooley fc Fender 1997 ; 
Foster et al. 1996[ ) show that GRS 1915+105 exhibits 
several giant flares per year, not all of which were ob- 
served with detailed campaigns (e.g., Fender et al. 1999a ). 
This yields an estimated average kinetic power, and, 
since almost all of the energy will initially be deposited 
in the form of CRs, an estimated cosmic ray power of 
ikin > lO^ergss" 1 for GRS 1915+105 alone. In fact, 
GRS 1915+105 seems to release an even higher power 



Kuulkers 2001), giving a very conservative lower limit on 
the Galactic rate of 1 yr" 1 . As with GRS 1915+105, we 
expect many giant flares to have gone unnoticed, and a 
more reasonable estimate of the event rate would be of 
the order of 10-100 Galactic events per year. 

Cyg X-3, which is believed to be relativistically 
beamed, shows radio peak luminosities up to 200 times 
stronger than GRS 1915+105 (|Fender fc Kuulkers 200l|), 



and often exhibits flaring activity( pgley et al. 2001 ) at or 
above the peak level of GRS 1915+105 on timcscales of 
^10 days. The other sources mentioned abov e are very 
similar to each other in peak radio luminosity ( Fender fc 



mated to exceed few x 10 3Y ergs s 1 ( Mirabel & Rodriguez 
1999| ) and even 3 x 10 38 ergs s" 1 flFender fc Pooley 2000| ) 



in the |form of microflarcs between major outbursts, esti- Kuulkers 2001), which we take as an indicator of kinetic 



Using the publicly available GBI monitoring data 



(tittp : //www . gb . nrao . edu/f gdocs/gbi/gbint .html), 

we estimate that GRS 1915+105 spends in excess of 60% 
of its time at flux levels significantly enhanced over the 
baseline flux (GBI monitoring data of Cyg X-3 show 
a similar rate), with about 2 major outbursts per year 
(see Fig. || and also Fender et al. 1999a; Foster et al 



199E). Assuming that the observed radio flux in flares 



is proportional to the amount of kinetic energy released 



in the flare, and using the observed 1997 flare (Fender 



power (most of these sources are not resolved and an es- 
timate of the equipartition energy of the jet is thus not 
possible). 

If indeed these sources operate on the same level as 
GRS 1915+105, the minimum kinetic luminosity of these 
seven sources together would be Lkin ^ 10 39 ergss _1 . 
Since this estimate is based on the minimum energy esti- 
mates of I/kin in GRS 1915+105, the true kinetic luminos- 
ity of these sources might well be much larger. 

Furthermore, there are many sources that are known 
to have been active at earlier epoch s [e.g., V404 Cyg (Han 
fc Hjellming 1992|) and Cir X-l, QHaynes et al. 1978|)], 



et al. 1999a) with a minimum kinetic energy estimate of 
2 x 10 44 ergs, we estimate that over the period covered 
by GBI monitoring, the mean kinetic luminosity of GRS 
1915+105 in flares is of order I/kin ~ 10 38 ergss _1 . If the 
baseline radio emission from GRS 1915+105 is also due 
to low level jet emission, this estimate increases by a 
factor ~ 1.4. 

The jets in SS433 are even more impressive: 
Reasonable estimates put the total, continuous kinetic 



which exhibited flux levels comparable to the above men- 
tioned sources. Many sources currently active might sim- 
ply not have been detected yet. Similarly, many more X- 
ray sources are observed to be consistently active at lower 
radio fluxes (e.g., Cyg X-l, or GX 339-4, |Fender 200i| ) 
than the brightest sources mentioned above. 

During the past few years, it has become clear that 
radio emission from Galactic X-ray sources is a very com- 
mon phenomenon. Radio emission is usually detected dur- 



Heinz & Sunyaev: A Narrow Component to the CR Spectrum 



1.00 F 






500 1000 
MJD - 50000 



1500 



Figure 6. Plot of the 2.25 GHz GBI monitoring data ( ittp : //www . gb . nrao . edu/f gdocs/gbi/gbint . html ) for GRS 
1915+105 over the time span from June 1994 to August 2000. Shown as a hatched area is the flux considered above 
the baseline flux, i.e., the flux considered to originate from flares, which we integrated to arrive at the estimate for the 
average kinetic power carried by the jet. The flare analyzed by Fender et al. (1999a ) is indicated by the mark "F99". 
The large gaps are due to gaps in the monitoring campaign and were not included in the procedure. 



ing state changes of the X-ray source (into or out of the 
low/hard state), including soft X-ray transients. While the 
powerful radio flares discussed above are associated with 



to extend only over the range in frequency observed in the 
radio. While an IR detection of the jet indicates a high en- 
ergy t ail of the spectrum flSams et al. 1996 ; Mirabel et al 



such transients, there arc many more X-ray sources which 1996; Eikenberry et al. 1998), a low energy component 



are active at lower radio flux levels (e.g., Cyg X-l, or GX 
339-4, [Fender 2001| ). 

These sources are observed to produce stationary, op- 
tically thick jet emission (as opposed to the already opti- 
cally thin emission detected in typical radio flares of tran- 
sient jets). It is not clear whether these jets are in fact 
relativistic and how much energy they carry. One might 
hope to estimate the kinetic power from the observed flux, 
scaling it to the peak flux observed in GRS 1915 as was 
done above for transient sources, but detailed kinematic 
modeling of the jet would be required to justify such a 
simple argument. In any event, because no complete sam- 
ple of such sources exists, it is impossible to estimate 
the total fraction of mechanical jet power contained in 
low power sources. Furthermore, these sources might have 
shown transient activity in the past as well, given that 
GRS 1915+105 also shows steady, quiescent radio emis- 
sion at comparable flux levels. 

Other sourc es like IE 1740-294 ( |Mirabel et al 1992|) 
GRS 1758-258 ( [Marti et al.^99^ ), and Cir-Xl ( [Stewart 
et al. |.995 ) show permanent extended structure resem- 
bling radio lobes in extragalactic radio sources, which are 
witness to past radio activity and must harbor a signifi- 
cant amount of CRs as well. We note here that estimating 
the total kinetic power from the presence or absence of ra- 
dio lobes in microquasar sources (indicating past activity) 
is severely hampered by the fact that Galactic radi o lobes 
are expected to have very low surface brightness flHeinj 
20020 

The estimate for the kinetic energy output from GRS 



1915+105 (Fender et al. 1999a), which we used as a tem- 
plate case to estimate the total Galactic energy in jets, is 
based on the assumption that the jet plasma is composed 
of a powerlaw of relativistic electrons and cold protons (for 
charge conservation). The electron spectrum was assumed 



(down to 7 ~ 1 — 10) has never been observed in any jet 
due to lack of viable emission mechanisms to reveal such 
a component. The possibility of detecting this component 
via inverse Compton scattering has been discussed, for ex- 



ample, by Ensslin & Sunyaev (2001). 

Finally, the physical structure of microquasar jets is 
still not known - they could be made of either discrete 
ejections or a continuous stream of matter. If the jet is not 
composed of discrete ejections, but instead is a continu- 
ous outflow with knots corresponding to internal shocks, 
Kaiser et al. (200"o| ) have shown that, in the case of GRS 
1915+105, the estimate of the total kinetic energy carried 
in the jet (and thus the total CR energy released in the 
working surface) is a factor of 10 higher than the above 
estimate (though the instantaneous kinetic power is re- 
duced), corrections for the low energy end of the particle 
distribution notwithstanding. 

All of this indicates that the lower limit of L^ n ~ 
few x 10 38 ergs s _1 is conservative, and it might be that 
the kinetic luminosity from microquasars is of the order of 
10% of the total Galactic cosmic ray power Lcr. Clearly, 
the uncertainty in frequency and power of radio flares 
in microquasars warrants continued monitoring of these 
sources to answer the question of how important energy 
input by these sources really is. 



3.2. The termination shock 

In the following, we will assume that a standard working 
surface as depicted in Fig. exists. The conversion of ki- 
netic to internal energy will then take place either in the 
forward or the reverse shock. 
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3.2.1. Discrete ejections 

If the jets are composed of discrete ejections, then each 
ejection will eventually convert its kinetic energy into 
some form of internal energy via interaction with the ISM, 
not unlike the external shock encountered in GRB after- 
glows. During the phase of relativistic propagation of the 
ejection, either the forward or the reverse shock (or both) 
must be relativistic. It is inside of the relativistic portion 
of the shock that most of the energy is dissipated. 

For a cold ballistically expanding ejection, most of the 
energy is dissipated in the forward shock. This is implied 
by the small opening angle 9 of the ejection: as long as 
9 r~j Cs/Tjctc <C 1/Tj e t (with Tjct being the bulk Lorentz 
factor of the ejection), the characteristic transverse size of 
the ejection R is always much smaller than the distance d 
over which it slows down, as seen in the frame of the blob: 
R ~ 9d <C d/T j e t (here, the factor of l/Tjet accounts for 
the Lorentz contraction in going to the frame of the ejec- 
tion). Thus, the deceleration time d/Tj ct c is much longer 
than the light crossing time of the ejection R/c, and the 
deceleration must occur in a sub-relativistic shock. This 
implies that the ejection is not heated to relativistic tem- 
peratures, while the forward shock must be relativistic, 
with shock velocity corresponding to rj et . 

The particles passing through this shock will have en- 
ergies of order rj et m p c 2 < rj e t,o m pC 2 , where 1^,0 is 
the initial Lorentz factor of the ejection (before interac- 
tion with the ISM). Because the ejection is slowed down, 
particles of a spread in energies are created in the forward 
shock, though the spectrum will show a sharp turnover or 
cutoff at energies e J> Tj eti o rn p c 2 (below this cutoff, it is 
plausible that the spectrum rises with dN/dj oc 7 -3 , see 
Appendix [X]). 

3.2.2. Continuous jets 

However, there are reasons to believe that at least in the 
well studied cases of GRS 1915+105 ( |Kaiser et al. 2000| , , 
though this question is up for debate) and in SS 433, as 
well as in the persistent radio outflows of low-hard state 
sources like GX 339-4 and Cyg X-l the jets are contin- 
uous streams of matter. In which part of the shock the 
dissipation occurs in a continuous jet depends on the jet 
thrust and the external density: The advance speed of the 
jet is roughly given by ram pressure balance. The transi- 
tion from forward to reverse shock dissipation then occurs 
when 



2L 



^-shock Py. C 3 



'16 ^0.1 n ISM 



(2) 



where L 38 . 5 = L iin /10 ds a ergs s - L kin /3 X 
10 38 ergs s _1 is the kinetic power, l 16 = ^ ct /10 16 cm 
is the jet length, j4 s hock the surface area of the shock, 
#o.i = 9/0.1 the jet opening angle in units of 0.1 radian or 
5.7o, and uism = Px/m p the ISM particle density in units 
of 1 cm~ 3 . Note that the working surface propagates, so 
he is a measured quantity, given by the age of the source, 



its power, and the ISM density. Strictly speaking, the jet 
thrust will have to be measured in the shock frame, and 
the relativistic shock jump conditions applied, but in or- 
der to estimate which of the two shocks is going to be 
relativistic, the above approximation is sufficient. 

The reverse shock is relativistic if the ratio on the left 
hand side of Eq. (j^) is smaller than unity, the forward 
shock is relativistic if the ratio on the left hand side is 
larger than unity. For Galactic sources both cases can oc- 
cur for appropriate external densities, depending on the 
length of the jet. 

3.3. Particle acceleration in relativistic shocks 
3.3.1. Nucleon acceleration 

For non-relativistic shocks, it has long been known that 
multiple shock crossings can accelerate particles to ultra- 
relativistic energies. The most commonly accepted scheme 
is first order Fermi acceleration (e.g., Krymskii 1977; 



Blandford & Ostriker 1978; Bell 197S; Blandford & Eichlci 



1987 ), which typically produces a powerlaw distribution 

/(7) = /o7" s (3) 

with canonical index s ~ 2. 

For relativistic shocks the situation is not as clear cut. 
Several attempts have been made to solve the problem 
of particle acceleration at relativistic shocks, mostly in 
the limit of test particle acceleration. In general, it is 
found that powerlaw distributions with somewhat steeper 
spectra than in non-relativistic shocks can be produced 



(Ellison & Double 2002; 


Achterberg et al. 2001 


Schneider 1987 


)■ 



It is certain, however, that acceleration of particles in 
the shock must take place: Particles crossing the shock 
are by nature already relativistic in the downstream rest 
frame, with a typical Lorentz factor of r rc i, the relative 
Lorentz factor between upstream and downstream frames. 
The particle distribution leaving the shock is thus strongly 
anisotropic, and essentially mono-energetic. The random- 
ization of this energy is then a question of the efficiency of 
the typical plasma processes often assumed to be present 
in populations of relativistic particles. 

The simplest assumption is that the particle momenta 
are simply isotropized behind the shock. The shock ac- 
celeration kernel is then a delta function and a cold up- 
stream plasma will be transformed into a narrow but rela- 
tivistic energy distribution, the width err of which should 
roughly be given by the Lorentz transformed width a of 
the upstream energy distribution, cjy ~ r rc i/3 a, where 
r rc i is the relative Lorentz factor between the upstream 
and downstream frames. The mean particle energy will be 
(7 m (?) ~ F ro i m c 2 ~ Tjct m c 2 . 

If scattering by downstream turbulence or particle in- 
teractions is stronger, the particles can be thermalized, in 
which case a Maxwell-Boltzmann distribution according 
to the relativistic Rankine-Hugoniot jump conditions will 
be established. The main observational difference between 
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these two cases will be the width of the energy distribution 
(see Fig. ||). 

If a significant fraction of the particles can perform 
multiple shock crossings (which again hinges on effective 
scattering to return downstream particles to the shock), 
we expect a powerlaw-type distribution to be established. 
It is reasonable to assume that the first time escape frac- 
tion from the downstream region (i.e., the probability that 
a particle which crossed the shock only once) is of order 
■Pesc <J 90% (e.g., Achterberg et al. 2001), which implies 



that most of the particles will only cross the shock once 
(note that the escape probability in non-relativistic shocks 
is generally very small for fast par ticles, tho ugh V esc is of 
order unity for thermal particles, Bell 197S ). The escape 
fraction in the upstream region is much smaller and gen- 
erally neglected in calculations. 

The particles re-crossing the shock will pick up another 



factor of order r 2 e j in energy gain ( Vietri 1995), which im- 



plies that a significant fraction of particles will be boosted 
to higher energies (7 > r 3 ol ). This fraction of the particles 
will contribute a significant amount of pressure to the post 
shock gas, which will modify the shock structure accord- 
ingly. Thus, the amount of energy accessible to the bulk 
of the particles which cross the shock only once is of the 



order [1 + r 2 cl (l/7? csc - 1)]" 



50% - 90%. 



The remaining fraction 1 — V csc of the particles will 
perform true Fermi acceleration. The low energy turnover 
of this distribution should be located roughly at 7 ~ r 3 cl 
( Victri 1995) and subsequent shock crossings will produce 
features at energies of r 2 ^ +1 (where i enumerates the num- 
ber of shock crossings). The similarity of this process to 
Compton upscattering was already mentioned in Sect. 2.3. 

The superposition of features from multiple shock 
crossing cycles will lead to the formation of a powerlaw 
at high energies, very similar to the powerlaw produced 
by optically thin inverse Compton scattering (where the 
Lorentz transformations of the photon distribution to and 
from the particle rest frame are replaced by Lorentz trans- 
formations to and from the upstream fluid rest frame, as- 
suming that scattering is strong enough to isotropize the 
particle distribution. The optical depth tic is replaced by 
the return probability 1 — Vesc-) The shock powerlaw s lope 
Sshock is determined by r ro i and P esc (e.g., Bell 1978): 



-^shock ^ 1 



In (1- Pose) 

In (r 2 ) 



(4) 



Thus, observationally, diffusive acceleration will lead 
to the presence of a second feature in the CR spectrum 
around ~ r 3 et . This situation is shown in Fig. where 
the powerlaw spectral component was plotted for Tj Ct = 5. 
The energy of this broad peak coincides nicely with the 
peak in the excess CR proton component required in the 



HEMN model ( ptrong et al. 2000| ) to explain the 7-ray 
EGRET excess above 1 GeV. 



We note that Achterberg et al. (2001) argue that in 
the absence of an efficient scattering mechanism, the av- 
erage energy gain per particle will be restricted to a factor 
of order unity (rather than r 2 e j) for higher order shock 



crossings (i > 2), and that the low energy turnover of 
the powerlaw component might thus be located at ~ T 2 cl 
(rather than r 3 cl ). The search for additional features in 
the CR spectrum at energies T 2 and T 3 might offer a po- 
tential way to test these predictions. 

It is also possible that a population of relativistic pro- 
tons already exists upstream of the shock. If the proton 
number density is equal to the electron number density, 
the bulk of the particles and of the energy will presum- 
ably be at the low energy end (7 ~ 1), otherwise the 
estimates for kinetic energy flux in the jet would have 
to be increased accordingly, increasing the impact on the 
Galactic CR spectrum as well (by a factor of 70, the low 
energy cutoff of the distribution) . 

This population will be shifted to higher energies in the 
shock, and possibly experience further Fermi type accel- 
eration. Thus, if a powerlaw of relativistic protons exists 
prior to the terminal shock with lower cutoff 70 and spec- 
tral index Sj et , the terminal shock should shift the lower 
cutoff roughly to rj Ct 7o, while the slope of the new power- 
law will be Min(sj 0t , s s hock)> where s s hock is the powerlaw 
slope produced in the relativistic terminal shock. 

While the particles produced in the shock must be 
relativistic at injection, the dynamical evolution of the 
shocked gas can reduce their energies significantly if adia- 
batic cooling is important before the particles can escape 
the shock region (radiative cooling of the nucleon distri- 
bution will be negligible) . The diffusion of particles out of 
radio lobes and hot spots is a highly uncertain process and 
has not been studied in the necessary detail to answer this 
important question. Rather than discussing it at length, 
which would by far exceed the scope of this paper, we de- 
cided to include a short discussion in Appendix |b], where 
we show that adiabatic losses do indeed pose a significant 
obstacle to particle escape (see also Fig. ||). 

3.3.2. Electron acceleration 

CR protons do not emit any measurable amount of syn- 
chrotron radiation. Inverse Compton scattering on protons 
is also very inefficient. However, if diffusive acceleration 
is operating efficiently in the working surface, then the 
powerlaw electrons produced there should be detectable 
through their radiative signature. Tentative evidence of 
this powerlaw electron component are the radio lobes in 
the microquasars IE 1740-294 and GRS 1758-258. 

The spectral index of the non-thermal emission from 
the lobes of IE 1740-294 is a ~ 0.7 - 0.9 flMirabel et aT 



1992), corresponding to an electron spectrum of slope of 
s ~ 2.4 — 2.8. If the protons are injected at the working 
surface with the same slope (this would require that the 
electron spectrum is unaffected by radiative cooling and 
that the powerlaw electrons are not just advected through 
the shock from upstream), then the observed CR slope 
near the earth should be steepened by As = 1/2 due to 
the energy dependence in the Galactic diffusion coefficient. 
This would yield a proton powerlaw slope of s ~ 2.9 — 3.3, 
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somewhat steeper than the canonical CR slope of s ~ 2.7, 
though not significantly. 

Part of the observed powerlaw electron distribution 
inside the tentative radio lobes of IE 1740-294 and GRS 
1758-258 could also have been advected from the jet. In 
fact, is is unclear whether electrons can be accelerated ef- 
ficiently in a shock if protons are present, since the shock 
thickness will be set by the proton Larmor radius, which 
will be much larger than the electron Larmor radius (e.g., 
Achterberg et al. 2001). In this case, the electrons will not 
experience a shock at all, more likely, they will be acceler- 
ated adiabatically to a narrow component with energies of 
order Tj ct m c c 2 . Such a component will not be detectable 
through synchrotron radiation. 

A detailed model of the energetic history of the par- 
ticles in the lobes of IE 1749-294 would be required to 
answer these questions. 

4. Observational consequences 

4.1. Predicted contribution to the Galactic CR 
spectrum 

Based on the picture laid out in Sect. |3|, we can try to 
predict what might be observed in the CR spectrum due 
to the presence of Galactic microquasars. 

First, it is obvious that close enough to a powerful 
relativistic jet source the locally observed CR spectrum 
will be completely dominated by the CRs produced in the 
terminal shock of the jet. However, it is clear that the 
powerlaw spectrum observed near earth is not dominated 
by a narrow component of microquasar origin - the current 
spectral limits rule out any contribution greater than a few 
percent. 

In a simple isotropic diffusion picture, the CR energy 
density in the environment of a continuously active source 
will fall roughly like the inverse distance to the source r -1 
[see Eq. (|§|)] for large r much larger than a particle mean 
free path, r ^> k/c ~ 0.2 pc -ff 2 pc /ri5, and smaller than 
the Galactic disk height, r <; #disk = 1 kpcTJkpc- 

Given the observed CR energy density of ~ 
10 -12 ergs cm -3 , we can estimate that the sphere of influ- 
ence of a given source, defined as the region inside which 
the source contributes more than 30% of the total mea- 
sured CR power (at which level it would enter the realm 
of detectability of by A MS 02) has a radius of order 



R 



30% 



1 kpc- 



^38.5 
«28.3 



(5) 



The time it takes the source to populate this sphere with 
CRs is roughly 



r 30% ~ 15 x 10 6 yrs 



T 2 



^28.3 



(6) 



Particles that reach the edge of the Galactic disk will 
be siphoned off into the Galactic halo. Thus the CR 
flux will diminish rapidly, roughly exponentially, beyond 
a source distance of order the disc thickness, r > H. It 



would take source activity longer than 10 s yrs to fill the 
Galactic halo and the disk. 

The well known microquasars mentioned above are all 
located much further from the solar system than this limit. 
However, if a source similar to, say, GRS 1915+105 had 
been active in the solar neighborhood (inside about 1 kpc) 
within the last ~ 10 7 yrs, our local CR flux should show 
a clear sign of the contribution from this source. 

In this context it is important to mention that GRO 
J1655-40, V4641 Sgr, Cyg X-3 (and also SS433) are known 
to be in high-mass X-ray binaries. Their lifetimes are 
therefore expected to be short. If such a relativistic jet 
black hole binary was located in the Orion nebula region 
within the past 10 6 yrs, we should be able to detect a 
strong signal in the low energy CR spectrum from this 
source alone. 

Far enough away from any single source, an observer 
will measure the time averaged contribution from all 
Galactic sources, washed out by CR diffusion (similar to 
the situation described in Strong & Moskalenko 2001). 



Since sources will likely follow a distribution of Lorentz 
factors of width Arj et , the observed signal will be smeared 
out over at least that width. Any intrinsic width of the 
produced CR spectrum will add to this effect, as well as 
broadening effects like solar modulation and scattering off 
of interstellar turbulence. 

In Fig. |?] we have plotted possible contributions to the 
CR proton spectrum from a single Galactic jet source. 
Depending on how much we have underestimated the 
power in Galactic jets and how much adiabatic losses of 
particles trapped in adiabatically expanding shock will 
suffer, we might over or underestimate the contribution. 
Taking the figure at face value, however, it seems likely 
that a contribution at the few percent level can be ex- 
pected in the energy region of a few GeV. 



4.2. Detectability of narrow features in the CR 
spectrum 

Detecting and positively identifying such a CR component 
will be a formidable challenge. The advent of the high 
sensitivity AMS 02 instrument and of the solar minimum 
might make it possible, however. Given the preliminary 
specifications of AMS 02, we can estimate the detectabil- 
ity of spectral features such as produced by microquasars. 
The rigidity resolution (i.e., energy resolution) of the in- 
strument is expected to be around 2% in the crucial range 
from 1 to 10 GeV, which will easily be sufficient to identify 
and resolve even the narrowest feature in Fig. ^. 

For an effective area of order 0.4 m 2 sr, the expected to- 
tal CR proton count rate by AMS 02 in the energy range 
from 1 to 10 GeV should be of the order of 10 3 s _x . At 2% 
energy resolution, this implies a detection rate of about 
2 x 10 8 yr -1 bin -1 , with a relative Poisson-noise level of 
order 10 . Calibration and other systematic errors will 
likely dominate the statistics, however, these numbers are 
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Figure 7. Toy model of the microquasar contribution to 
the CR spectrum, for a single microquasar situated in a 
low mass X-ray binary, active for t > 1.5 x 10 7 yrs on 
the level of 3 x 10 38 ergs s" 1 (similar to GRS 1915+105), 
and at a distance of 1 kpc. For simplicity, we assumed the 
source was operating with uniform bulk Lorentz factor 
Tjct = 5 (top) and I^t = 2.5 (bottom). The curves are 
normalized relative to the measured differential Galactic 
CR background spectrum (thick grey dashed line). Shown 
are the same curves as in Fig. ^ narrow feature for 
upstream temperature of T ~ 7 x 10 10 K {dotted line) 
and for T ~ 7 x 10 s K [solid line), Maxwellian with 
kT ~ rjetTO p c 2 /3 (dashed line), multiply scattered com- 
ponent for efficient pitch angle scattering (dash-dotted 
line) and for inefficient pitch angle scattering (dash-triple- 
dotted line), and the relative contribution of heavy ele- 
ments for metallicity 10 times the solar value (long-dashed 
grey curve), as seen in GRO J1655-40 and V4641 Sgr. 
Each spectral component has been steepened by E~ x / 2 to 
account for the energy dependence of the diffusion coeffi- 
cient. 

encouraging, and we expect that a source at the few- 
percent level will be detectable with AMS 02. 

The heavy element sensitivity of AMS 02 will share 
similar characteristics: for the same energy resolution and 
effective area, the detection rates of Carbon and Iron, 
for example, should be of order 4 x 10 s _1 bin -1 and 
4 x 10 4 s _1 bin -1 respectively. Aside from AMS 02, signa- 
tures might be detected by other instruments, and even ex- 
isting data sets might contain signals. Identification would 
require scanning these data with high spectral resolution. 
Note that the effects of solar modulation will broaden 
any narrow spectral component significantly. Results by 



~ 5 GeV will be broadened by ~ 1 GeV, (less at higher 
energies) though this effect will be reduced at solar mini- 
mum. 

4.3. Gamma ray emission from pion decay 

As the CRs produced in microquasars travel traverse the 
Galaxy, they will encounter the cold ISM. The interaction 
of a CR proton (by far the most abundant and thus most 
energetic component of the CR spectrum) with a cold ISM 
proton can lead to secondary particle production and to 
the emission of gamma rays via several channels, the most 
important of which is 7r° decay. 

Using the toy model presented in Fig. 0, we can es- 
timate how much gamma ray flux can be expected from 
the CR halo of a powerful microquasar and compare it 
to the background flux from the Galaxy. We assume that 
the CRs diffuse away from the source until they reach the 
Galactic halo, approximated as a zero pressure boundary 
condition at radius R ~ i?disk = IkpciJkps (assuming 
spherical symmetry for simplicity) . The result is shown in 
Fig. | 

Note that the gamma ray signal even for a source of 
average kinetic power of Lkin = 3 x 10 38 ergss~ 1 is small 
compared to the background signal coming from the same 
solid angle (nR 2 ). However, because the CR density in- 
creases towards the center of the source, higher spatial 
resolution can improve the signal-to- noise ratio somewhat. 
For a spherically symmetric cloud of CRs with luminosity 
L and vanishing pressure at the boundary R ~ -f/disk, the 
density follows 



L 



= 6 x 10 _11 cm- 3 



J d"f'"f'mc 2 f(y) 4irn 
f(l) 



1 1 

r ~ R 

i38.5 T 15 



(7) 



/ dyymc 2 f(y) h 



2 

kpc 



* 1 

Rkp<* V T 



Obviously, then, the gamma ray surface brightness of 
such a source is centrally peaked and thus the best ob- 
serving strategy is to try to resolve the source. The photon 
flux F v from within an angular distance < r / D from the 
source position (here, D is the physical source distance, 
while the on-sky angular distance from the source posi- 
tion) is then proportional to 



F„ 



<-) 



AttD 2 3 k 
x |l + 3r 2 

fl 



J dyymc 2 f(y) 



lnfl + Vl^ 



In (r) 



(1 + r 2 ) Vl ~r 2 } 



(8) 



whereas the galactic background flux F v ^ back for the same 
solid angle in the small angle approximation is simply pro- 
portional to (r/D) 2 . Thus, the source contribution relative 
to the background, measured as F Vt a omce/Fv back) rises to- 
wards lower r, which we have plotted in the insert in Fig. ^| 



Labrador & Mewaldt 1997 demonstrate that a line at 
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8/D ln [arcmin] 

"to 




them from traveling further than about 250 pc from the 
source. 

However, they could act as a significant ionization 
source for the surrounding medium: The ionization loss 
timescalc for a particle with energy E = 30 MeV E30 is of 
order flGinzburg 1979| ) 



1000.0 



Figure 8. Toy model for the gamma ray signature pro- 
duced in a microquasar CR halo via pion decay (includ- 
ing 7T° decay and bremsstrahlung from secondary electrons 
and positrons). Curves were calculated assuming a source 
active for over r ^> 15 x 10 6 yrs with an average power of 
-kkin = 3 x 10 38 ergs s _1 , at a distance of 10 kpc, and for an 
ISM particle density of 1 cm" 3 , assuming CRs are lost once 
they have reached the edge of the Galactic disk at about 
1 kpc distance from the source. Labes according to Fig. ||. 
dotted line: Gamma ray signature from narrow component 
of the microquasar halo; dashed line: from thermalized 
component; dash-dotted line: from " powerlaw type compo- 
nent" in the case of efficient scattering; dash-triple- dotted 
line: from "powerlaw type component" in the case of inef- 
ficient scattering. Thick grey line: EGRET diffuse gamma 
ray background at the position of GRS 1915+105 (Hunter 
et al. 1997), assuming the same solid angle as subtended 
by the source. Thick long dashed line: background gamma 
ray emission over the same solid angle modeled assuming 
the proton CR spectrum measured near earth and an av- 
erage ISM density of 1 cm -3 . Hatched line: GLA ST sensi- 
tivity. Models were computed using the GALPROP routines 
by Moskalenko & Strong (1998). Insert: Angular depen- 
dence of the source contribution to the gamma ray flux 
F v , relative to the background flux B v at a photon energy 
of 1 GeV [see Eq. (|j)]. The right Y-axis shows the source 
distance at which the GLAST sensitivity is reached for the 
value of F v jB v shown in the curve, the top X-axis shows 
the angle 8 = r/D in units of arcmin/Dio where -D10 is 
the source distance in units of 10 kpc. 



4.3.1. Sub-cosmic rays from sub-relativistic sources 
like SS433 

Objects like SS433 might be important sources of sub- 
cosmic rays in the Galaxy: The cold thermal ions carried 
in the jet at 0.26 c will be accelerated to energies of order 
E- lon = x 30 MeV/nucleon. These particles will suffer 
severe ionization and Coulomb losses which will prevent 



E. 



3/2 



2 x 10 6 yrs— 22-, 

nisM 



(9) 



much shorter than the expected lifetime of the particles 
inside the Galactic disk, i.e., the particles will deposit their 
energy in the vicinity of the source, especially if the source 
is located within a molecular cloud. 

Furthermore, the excitation of nuclear 7-ray emission 
lines by interaction of these sub-cosmic rays with inter- 
stellar heavy ions of C, O, Fe, and other elements might 
be detectable by INTEGRAL. 

4.4. Cold electron-positron jets 

If the jet consists chiefly of relatively cold electron- 
positron plasma, and if dissipation occurs mostly in the 
reverse shock, then the jet terminus will produce relativis- 
tic electrons and positrons with energies of the order of 
Tjet m e c 2 ~ 2.5 MeV, which will then begin to diffuse into 
the ISM. Such positrons and electrons could produce addi- 
tional bremsstrahlung radiation at energies of a few hun- 
dreds of keVs up to 2.5 MeV. Much like mildly relativistic 
protons, these electrons will contribute to the heating of 
the ISM due to the ionization losses, but much more im- 
portantly for future observations might be the electron- 
positron annihilation line at 511 keV. 

For an integrated mechanical luminosity of L^in = 
3 x 10 38 ergs s" 1 L38.5 of the entire ensemble of Galactic 
relativistic jets, the flux of positrons carried by jets is 



v kin 



2 (7) Fjet m e c 2 



4 x 10 43 s _1 



^38.5 

r 5 ( 7 )' 



(10) 



where (7} is the mean random Lorentz factor of the elec- 
trons/positrons in the jet plasma. Note that this esti- 
mate will be independent of whether dissipation occurs 
mostly in the forward or reverse shock. Assuming that 
these positrons can diffuse out of the radio plasma they 
are originally confined in, the estimated luminosity of the 
annihilation line will be of the order of L a nn ~ £jet/(7}- 

This is actually comparable to the total amount 
of positrons annihilating in the Galaxy according to 
the observations of the e + /e" annihilation line from 
OSSE/GRO, N e+Gal ~3x lO^s" 1 (|Purcell et al. 1997|). 



If Galactic jets are in fact composed of electron-positron 
plasma, this measurement immediately implies one of the 
following conclusions: a) either the mechanical luminosity 
of these jets is not far above our relatively conservative 
estimate of 3 x 10 38 ergss _1 , or b) the pair plasma is not 
cold, i.e., (7) ^> 1, or c) diffusion of particles across the 
magnetic boundary of the remnant jet plasma is very inef- 
ficient, in which case many Galactic "radio relics" should 
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exist, not unlike in the case of radio relics from radio 



loud AGNs in the intracluster medium (e.g., Ensslin ct al 
19980 

The Integral SPI spectrometer and the IBIS imager 
would be able to measure the increase in the annihila- 
tion line flux towards microquasars located away from 
the Galactic center (where the background is highest) like 
GRS1915+105, and to measure the line width if it could 
be detected. These measurements could be very helpful 
in constraining the particle content of relativistic Galactic 
jets. 

5. Conclusion 

We have argued that Galactic microquasars produce cos- 
mic rays with energies of a few GeV/nucleon. Under 
preferable conditions (i.e., if particles in the jet are 
cold and do not suffer significant adiabatic losses be- 
fore escaping the acceleration region) this component 
should be detectable as a narrow peak with energy ~ 
Tjet m p c 2 /nucleuon in the Galactic CR spectrum. The su- 
perposition of several such peaks would then reflect the 
distribution of Lorentz factors of Galactic jet sources. 

In addition, diffusive acceleration of particles might 
produce a powerlaw distribution of particles with a low 
energy turnover at energies around T? ot m p c 2 to r? et m p c 2 , 
visible as an edge-like feature in the CR spectrum. 

The locally measured contribution to the CR spectrum 
will be strongly dominated by sources operating close by 
(within a distance of about one Galactic disk height and 
within the past 10 7 yrs), since at distances much larger 
than that the contribution from a given microquasar falls 
off exponentially. We estimate the global energy content 
in the CR component accelerated in Galactic relativistic 
jets to be at the 0.1% to 10% level of the total Galactic 
CR luminosity. 

This CR contribution from Galactic relativistic jet 
sources might be strongly overabundant in heavy ele- 
ments, reflecting the composition of the accretion disk 
where the jets are launched. Thus, it is possible that the 
chemical abundance measured in the GeV region (where 
we expect the contribution from jets to show the strongest 
effect) will differ slightly from the chemical composition at 
higher energies. 

While signatures of the CR component from micro- 
quasars might already be buried in existing data, the up- 
coming solar minimum and the launch of AMS 02 will 
offer ideal conditions to search for this component and to 
put constraints on the microquasar activity in the nearby 
universe. 

We suggested that the absence of any observable traces 
of microquasar component in the cosmic ray proton spec- 
trum could be used to argue in favor of electron-positron 
jets. For this case, we showed that existing OSSE/GRO 
observations of the Galactic electron-positron annihila- 
tion rate can be used to limit the power in cold electron- 
positron jets inside the Galaxy to L± < 3 x 10 38 ergs s _1 . 
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Appendix A: The particle spectrum produced in 
decelerating ejections 

It is straight forward to calculate the particle distribution 
produced by single pass shock acceleration in a decelerat- 
ing ejection in the ultra-relativistic limit, assuming that 
one has knowledge of the single pass shock acceleration 
kernel at a given shock velocity r ro i. 

Take an ejection of initial mass M and Lorentz factor 
r , which is sweeping up and shocking external matter. 
The total energy of the ejection and the swept up matter 



E ~ Yc z \ Mn 



TV 



dN'T(N')m v 



(A.l) 



where the integration is over the number TV of swept up 
particles with mass m p at a given ejection Lorentz factor 
r, and the extra factor of T inside the integral takes ac- 
count of the shock acceleration of the swept up particles. 

S ince E is conserved, we can take the derivative of 
Eq. (|A.l|) with respect to T, and arrive at 



(J N dN'T(N')m T 

dT ~ 



Mn 



E 



T 2 C 2 



(A.2) 



for 1 <C T < Tq. For r ~ 1, the non-relativistic corrections 
reduce the amount of energy available, leading to a low 
energy cutoff at T ~ 1 and a high energy cutoff at r . 

To arrive at the observed particle distribution dN/dj, 
this must be convolved with the single shock acceleration 
kernel 7(r), however, for a narrow kernel, such as assumed 
in this paper, the powerlaw approximation seems suffi- 
cient: dN/dj cx 7~ 3 . 

Appendix B: Cosmic ray diffusion and adiabatic 
cooling 

An important question is whether the particles produced 



in the shock discussed in Sect. 3.3.1 can indeed diffuse out 



of the shock region, in which case they will freely escape 
and propagate through the Galaxy essentially with the 
energy obtained in the shock, or whether they are trapped 
inside the shocked gas until it expands adiabatically after 
the shock has passed and activity has ceased. In the latter 
case, the particles will lose a significant amount of energy 
to adiabatic expansion. 
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Following the discussion in Sect. 3.3.1, wc distinguish 
two cases: dissipation in the forward and in the reverse 
shock. 

The escape time of the particles out of the shock can 
be estimated as 



hock 



(B.l) 



where -Rshock is the typical size scale of the shock and k 
the relevant diffusion coefficient. 

While the shocked ISM must still be magnetically con- 
nected with the unshocked ISM, the jet plasma will be sit- 
uated on field lines advected out from the central engine, 
which are likely not connected with the ISM. In the for- 
ward shock, the relevant diffusion coefficient should then 
be taken as ^forward ~ ^ih the diffusion coefficient parallel 
to the mean magnetic field, while for the reverse shock one 
has to consider diffusion across the magnetic boundary of 
the contact discontinuity between shocked jet plasma and 
shocked ISM in addition to diffusion to the contact dis- 
continuity and away from it. 

A lower limit on the diffusion time out of the shock is 
thus given by the value for the forward shock, since the 
particles which have diffused out of the reverse shock must, 
in addition, also propagate through the forward shock. 

B.l. Forward shock: 

Using the simple approx imate expression for ku (e.g., 
Kennel fc Pctschck 1966 ), and assuming a Kolmogorov 
turbulence spectrum for the magnetic field originating 
on scales of order the shock size i? s hock = ^jct & = 



10 15 cmZi 6 $o.i an d containing a fraction e = B 2 nrh /B 2 ot 



of the total magnetic energy, the parallel diffusion coeffi- 



cient for a particle with energy 7 m c 
written as 
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Here, Scatter is the scattering frequency for particles 
with gyro radius r& off magnetic turbulence with energy 
density kB^/8n and wave numbers of order k ~ t"q . 
Scattering can also be induced by collective resonant in- 
teractions of particles with the field, in which case the 
parameter e denotes the efficiency of this process. 

Writing the shock area as A s h oc k = ^shock = n x 
10 30 cm 2 #0.1 '16 g lVes an approximate hot spot pressure 

Pshock of 



Pshock 



kin 



^■shock C 



= 3xl(T 3 ergs cm" 3 L 38 .5 lw 6»^f .(B.3) 



If the magnetic field is at a fraction tp of the equipar- 
tition field, ip = B/B cq , we have B = ip B cq ~ 

0.5GipLlil/h 6 e .i. 



The comoving (i.e., measured in the frame of the 
shocked plasma) limit to the proton escape time t!, is then 
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If we write the Lorentz factor of the shocked ISM gas as 
r', then time dilation gives t\\ — V r!, in the observer's 
frame. 

If the turbulent velocity inside the region of interest 
is comparable to the expansion velocity, and if large scale 
turbulence is present (which was the underlying assump- 
tion in out estimate of e above), then turbulent trans- 
port could aid particle escape: In a simple mixing length 
approach, the diffusion coefficient can be approximated 
by Kturb ~ 1/3 Wturb ^turb, where ^turb is the characteristic 
turbulent velocity and / t urb the scale length of the largest 
scale turbulence. The turbulent transport time is then 



Tturb 



3R< 



hock 



rJi ^shock 

10 s 



0.1 K 



hock 



. , . , (B.5) 

Vturb H V t iturb 

If the shock region is expanding roughly ballistically 
as it propagates, the adiabatic loss time scale measured in 
the observer's frame is 

1=4^ = ^-3x10^6 (i 
7 P i?shock \l. 



Tad 



(B.6) 



P -Rshock 

If the forward shock is relativistic, we need to take 
time dilation into account when comparing r a d and Tji . It 
follows that the adiabatic loss time scale r ac i is longer than 
the escape time only if 



7-I/6 n 'jot t^2/3 _ A -,n-3 

e 99X3^5^0.1 ^7 I5 < 4 x 10 7 

Since both e and ip might be significantly smaller than 
10 -5 , it is not implausible that tu ;> r a d. In the case of 
strong turbulent transport, this condition simplifies to 



(B.7) 
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If the external pressure is parameterized as pism = 
10 _11 ergs cm" 3 p_n, the maximum reduction in particle 
energy possible through adiabatic losses is roughly 
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2/5 



where the left hand side of the inequality corresponds to a 
relativistic equation of state in the hot spot and the right 
hand side to a non-relativistic equation of state. This, of 
course, means that a large fraction of the energy injected 
originally into CRs could be lost adiabatically, reducing 
the amount of energy available (estimated in Sect. 2.1) by 
up to a factor of ~ 1000. 
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B.2. Reverse shock: 

In jets where dissipation occurs mainly in the reverse 
shock, the jet geometry is similar to AGN jets (i.e., the 
jets arc effectively reflected by the ISM, thus inflating a 
cocoon with spent jet fuel). 

Particles accelerated in the shock will eventually be 
advected out of the shock region and into the cocoon (see 
Fig. ||). The timescale for this process is 



Tadvcct 



-Rshock 
^advect 



3 x 10 4 sZi 6 6»o.i 



(— r 

\ ^advect / 



(B.10) 



where w a dvect <^ c/3, since the limiting downstream veloc- 
ity in a strong relativistic shock is c/3. 

The diffusion time towards the contact discontinuity 
is given by Eq. (BT). In order to enter the forward shock, 
they have to propagate across the magnetic boundary at 
the contact discontinuity, which introduces an additional 
cross-field diffusion term. The contact discontinuity will 
have a typical thickness of the order of the Larmor radius 
tq of the particles, thus the particles will have to traverse 
a region of size tq perpendicular to the field in order to 
cross the contact discontinuity (this approximation is valid 
as long as the parallel diffusion time over one coherence 
length of the field is longer than the perpendicular diffu- 
sion time across one Larmor radius, otherwise the perpen- 
dicular diffusion time across one coherence length should 
be used). The lower limit to the diffusion time across the 

r 2 

field is then ti > — . 

The perpendicular diffusion co efficient kj_ c an be ap- 
proximated as K|| (kB 2 /B 2 ) 2 (e.g., Parker 1965| ), i.e., 
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The total diffusion time out of the reverse shock should 
then be of the order of r ~ rii + tj_, which has a minimum 

when f e ip 1 / 3 L^i/T^) ~3x 10~ 3 , with r min - 2 x 

V / mill 

10 5 Zi6 6'o.i. Thus, even in the optimal case, the diffusion 
time out of the reverse shock is longer than the advection 
time unless w a dvoct ^ 0.1 c. 

Thus, it is rather likely that the bulk of the par- 
ticles are advected out of the shock and into the co- 
coon. The pressure driven expansion of a cocoon can be 
approximated by a simple spherically symmetric model: 
Dimensional analysis sugges ts that the size of the cocoon 
R c follows a simple scaling ( Castor et al. 1975] ) 

/ T t 3 \ 1/5 

R c ~( — ) . (B.12) 

VPISM/ 

Using this scaling to obtain order of magnitude esti- 
mates of the conditions inside the cocoon, we can write 
the cocoon pressure p c as 

1 



Pc 



■ PlSM-^c 



2 x lO-^rgscm-^-g 4 / 3 ^^, 



(B.13) 



where R\§ = i? c /10 16 cm ~ lie is the cocoon radius. The 
magnetic field is 



Bc~3x 10- 2 G tp i? 16 2/3 L\il 



(B.14) 



Comparison of Eq. ( B.13 ) with Eq. (BJ3) shows that a 
(relativistic) particle advected out of the shock into the 
cocoon will suffer adiabatic losses of the order 
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which is only weakly dependent on the source parameters. 

Once the particles are inside the cocoon, the escape 
time is again given by r ~ T|| + t± with the expressions 
for T|| and t± from Eq. (B.4) and Eq. (B.ll), though with 
different values: 



T|| - 7 x 10 8 si?i 6 
and 

t± - 4.5xl0 3 si?ie 
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The adiabatic loss timescale is simply 
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Since this is proportional to R 5 ^ 3 , while Tj| oc i? 10 / 9 and 
t_\_ oc i? 8 / 9 , the particles are more likely to escape out of 
large, old cocoons than out of small, young ones. 

B.3. Summary 

The diffusion of CRs out of radio lobes is clearly an impor- 
tant problem not only in the context of CR emission from 
microquasars, but also for extragalactic radio sources and 
their contribution to the CR flux in clusters of galaxies. 
Since the current level of understanding is still very rudi- 
mentary, further study of this aspect is necessary. Given 
the large uncertainties in these estimates, and given that 
we tried to provide conservative estimates whenever pos- 
sible, we feel that there is a good chance that the mecha- 
nisms outlined in this paper will work under realistic cir- 
cumstances and that a measurable CR contribution from 
microquasars can be expected. 
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